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Quantum computation with a complete graph of superconducting qubits has been recently pro¬ 
posed, and applications to amplitude amplification, phase estimation, and the simulation of realistic 
atomic collisions given [Phys. Rev. A 91,062309(2015)]. This single-excitation subspace (SES) ap¬ 
proach does not require error correction and is practical now. Previously it was shown how to 
implement symmetric nxn unitaries in a single step, but not general unitaries. Here we show that 
any element in the unitary group U(n) can be executed in no more than three steps, for any n. This 
enables the implementation of highly complex operations in constant time, and in some cases even 
allows for the compilation of an entire algorithm down to only three operations. Using this protocol 
we show how to prepare any pure state of an SES chip in three steps, and also how to compute, for 
a given SES state p, the expectation value Tr(pO) of any nxn Hermitian observable O in a constant 
number of steps. 

PACS numbers: 03.67.Lx, 85.25.Cp 


I. PRETHRESHOLD QUANTUM 
COMPUTATION 

There is currently great interest in the develop¬ 
ment of special-purpose quantum computing devices and 
methodologies that do not require full error correction 
and which are practical now. For example, D-Wave Sys¬ 
tems produces commercial quantum annealers based on 
superconducting circuits that solve an important class 
of binary optimization problems [Tj- However it is not 
known whether the D-Wave annealers can outperform 
conventional classical supercomputers [1, Q. An opti¬ 
cal approach Q that solves an arguably less important 
problem—sampling from the distribution of bosons scat¬ 
tered by a unitary network—but which is likely capa- 
ble of quantum speedup has also been investigated 
0. An approach called the single-excitation-subspace 
(SES) method, also based on supercondonducting cir¬ 
cuits, has been proposed Q. Here computations are per¬ 
formed in the n-dimensional SES of a complete graph of 
n qubits. We call these examples prethreshold, referring 
to the threshold theorem of fault-tolerant quantum com¬ 
putation, because they do not require exceeding fidelity 
and qubit-number thresholds before being applicable. 

A quantum computer chip implementing the SES 
method consists of a fully connected array of supercon¬ 
ducting qubits with tunable frequencies and tunable pair¬ 
wise couplings; an abstract representation is given 

in Fig. [I] It works by operating in a subspace of the full 
2"-dimensional Hilbert space where the Hamiltonian can 
be directly programmed. This programmability elimi¬ 
nates the need to decompose operations into elementary 
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FIG. 1. (Color online) Complete graph with n = 16. The 
vertices (open circles) are qubits and the edges (colored lines) 
are tunable couplers. 


one- and two-qubit gates, enabling larger computations 
to be performed within the available coherence time. The 
price for this high degree of controllability is that the 
approach is not scalable. However, a technically unscal¬ 
able quantum computer is still useful for prethreshold 
quantum computation and might even be able to achieve 
speedup relative to a classical supercomputer for certain 
tasks. The SES approach trades physical qubits and high 
connectivity for, in effect, longer coherence. This is a 
sensible trade for quantum computing architectures such 
as superconducting circuits, whose largest prethreshold 
problem sizes are limited by coherence time, not by the 
difficulty of introducing additional qubits. A realistic 
chip layout that provides space for the coupler circuits 
and avoids the crossovers of Fig. [T] is shown in Fig. [5] 
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FIG. 2. (Color online) Possible layout for the 16-qubit chip. 


II. CONTENT OF THIS PAPER 

A significant restriction of the SES method presented 
in Ref. Q is that the Hamiltonian programmed into the 
hardware is real and symmetric, whereas the most gen¬ 
eral Hamiltonian is complex Hermitian. If a target oper¬ 
ation has the form where A is a known real sym¬ 

metric generator matrix, then the unitary can be imple¬ 
mented in one step. This is the case when the unitary 
is symmetric {U = U^) and is reviewed in Sec. IHI A1 In 
that section we also provide an improved procedure for 
constructing the time-optimal SES Hamiltonian T-L cor¬ 
responding to a given generator A. 

However, a general element of the unitary group U(n) 
has the form with M complex Hermitian. This 

is the nonsymmetric unitary case (U ^ U^) discussed in 
sec.nmi We show there that any nonsymmetric element 
U G U(n) can be implemented in three steps, for any n. 

Applications of these techniques are given in Sec. IIVI 
In Sec. IIV Al we show how to simulate time-independent 
but otherwise arbitrary nxn complex Hamiltonians with 
an SES chip in three steps. In Sec. IIV HI we show how to 
prepare pure but otherwise arbitrary SES states in three 
steps. And in Sec. IIV Cl we explain how to compute ex¬ 
pectation values of arbitrary nxn Hermitian observables. 


III. SES IMPLEMENTATION OF UNITARY 
OPERATORS 

A. Single-step implementation of symmetric 
unitaries 


The basic single-step operation in SES quantum com¬ 
puting is the implementation of symmetric unitaries of 


the form U = with A real and symmetric Q. 

Therefore, a standard task in SES algorithm design 
and implementation is the construction of an optimal 
protocol—an SES Hamiltonian 'H and evolution time 
tqc —to implement that unitary. We assume here that the 
generator matrix A is known; if it is not then the clas¬ 
sical overhead for obtaining A from U must be included 
in the quantum runtime. (We also note that the genera¬ 
tor A = ilogU is not unique.) The optimal protocol for 
implementing a symmetric unitary depends on the func¬ 
tionality assumed of the chip, especially of the tunable 
coupler circuits. Here we assume that the experimentally 
controlled SES Hamiltonian can be written, apart from 
an additive constant, as 

H = ffmaxAT with - 1 < Ku> < 1, (1) 

which we call the standard form. In this case we are 
assuming that the couplings can be tuned continuously 
between —(?max and gmax, and that the qubit frequencies 
can be varied within a window of width 2g„iax about some 
parking frequency. Because we are free to change the 
overall phase of an SES state, we write the symmetric 
unitary as 

U = (2) 


where I is the nxn identity matrix, and then ignore the 
global phase The value of c is chosen to minimize 
the evolution time tqc, which is proportional to the angle 

9a = niax- c5ii'\. (3) 

ii' 


The K matrix in ([T]) is then given by 


K = 


A-cl 
9a 


and the evolution time is 

_ h9A 

^qc — 

^max 


( 4 ) 

( 5 ) 


Note that 9a is not bounded by 27r and can become arbi¬ 
trarily large. The global phase angle that minimizes 9a 
is 


mini Ail + max* A^ 

c=-^-, (6) 

which is proved below. Although we have assumed that 
the SES Hamiltonian % — gmax AT is abruptly switched 
on for a time tqc before being abruptly switched off— 
which is the fastest protocol—any SES Hamiltonian of 
the form Ti, = g{t)K such that j{g/h)dt = 9a niay be 
used instead. 

To minimize over c we consider two cases: In the 
first case max^i' \ Aiii\ occurs for an ojj-diagonal element 
of A, in which case the minimum value of 9a is indepen¬ 
dent of c (because c only affects the diagonal elements of 
the shifted matrix A — cl). Therefore we only need to 
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consider the second case where maxii/ \Aiii\ occurs for a 
diagonal element. The diagonal elements consist of points 

X S {^11, A 22 , ■ ■ ■ , Ann} (7) 

on the real number line, bounded between min^ An and 
maxiAii. Placing c at the midpoint of the smallest re¬ 
gion containing all the points in Q minimizes the largest 
distance \Aii — c\. 


B. Three-step implementation of nonsymmetric 
unitaries: ABA decomposition 

Our protocol relies on the matrix decomposition 

U = 0,e-^^0j, (8) 

where D is a real diagonal matrix and the Oi G 0{n) 
are real orthogonal matrices. This identity follows from 
the KAK decomposition of the Lie group U(n) [^. To 
obtain the Oi and D from U, we first compute 

X=UU^ = Oie-‘^^°Oj, (9) 

which is both symmetric and unitary. The real and imag¬ 
inary parts of X &re also separately symmetric. Then the 
unitarity condition 

(Rey — ilmy) (Rex-b limy) = J (10) 

shows that Re y and Im y commute and can be simulta¬ 
neously diagonalized. Oi is determined by a Schur de¬ 
composition of Re y, which always produces a real Oi 
(unlike the decomposition of y itself). Then and 

O 2 are obtained from Oj yOi and respectively. 

The three-step implementation for a nonsymmetric 
U G U(u) follows from the identity 

U = (11) 


on a laptop computer 0. The quantum runtime to 
implement a nonsymmetric unitary is 


t 


qc 


h{29A + Ob) 

5max 


(16) 


with 9 defined in ©• The generator matrices A and B 
in dm) are not unique. 

The ABA decomposition allows for the possibility of 
implementing highly complex operations in three steps. 
But this does not imply that an entire algorithm, com¬ 
piled into a single unitary, can be implemented in con¬ 
stant time, because the compiled unitary might not be 
known a priori, and there is classical overhead (1151) for 
computing A and B. More importantly, evaluating A 
and B for an entire algorithm would presumably be pro¬ 
hibitive when one is attempting to outperform classical 
computers. Furthermore, algorithms might include mea¬ 
surement steps that cannot be postponed to the end. 


IV. APPLICATIONS 


A. Hamiltonian simulation 

A useful application of (fTTI) is to 17 = where 

H is a given complex Hamiltonian. In this case we have 

e ^ = e e e , (17) 

with A and B given by (USD and da, where A is a di¬ 
agonal matrix containing t/h times the spectrum of H. 
This enables the fast simulation of any time-independent 
Hamiltonian with an SES chip [HI- 


B. SES pure state preparation in 3 steps 

In some cases it is possible to compile an entire algo¬ 
rithm down to only three steps. As an example we give 
an algorithm for preparing any (normalized) pure SES 
state of the form 


which we call the ABA decomposition. Here A and B are 
real symmetric nxn matrices. To derive (ED we express 
the target unitary in the spectral form U = Ve~'^^V\ 
where V is complex unitary and A is real and diagonal. 
Decomposing V using ([5D we have 

U = Oie-^^Oj e-*^02 , 

= (OiOj) (OiOj f ^ ^^2) 

which leads to (ED with generators 

A = OiDOj, (13) 

B = 0i0JA020J, (14) 

which are both real and symmetric. The classical runtime 
to obtain A and B is about 


IV') Oil*), ai = |ai|e*®% 0 < < 27r. (18) 

i=l 

Here |i) = |0 • • • 1^ • • • 0) is the *th SES basis state of the 
n-qubit processor. We proceed by giving a protocol with 
linear depth that is practical for small n, which is then 
subsequently compiled down to three steps. 

We start with the basis state |1), which is prepared 
from the system ground state |00 • • • 0) by a microwave 
pulse, and then apply the standard-form SES Hamilto¬ 
nian T-L = ffmaxATstar for a time tqc = TTh!y/ng^i^y^, with 



/I 

1 

2 

1 

2 



1 

! 

2 

0 

0 ••• 

0 

-^star — 

0 

0 ••• 

0 


VI 

0 

0 ••• 

V 


1.4xn^-^/xs 


(15) 
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FIG. 3. (color online) Occupation probabilities for a uniform 
weight state. Phases of the probability amplitudes at are not 
represented in this figure. 


FIG. 4. (color online) Occupation probabilities for a typical 
target state. In this example the components of maximum 
and minimum weight have indices imax = 2 and imin = 3. 


the adjacency matrix for a star graph with qubit 1 at 
the center (see Sec. IIIA of Ref. [3]). This produces the 
uniform state 


M = 0{n) steps a uniform weight state is created. The 
final operation Wdiag shifts the phases of the uniform 
weight state to that of (|20)) . The first step is: 


if)^ 


|1) + |2) H + |n) 


( 20 ) 


apart from a phase. 

If the occupation probabilities in the target state m 
are uniform, 


1. Find the components Imin and Imax with the small¬ 
est and largest weights, respectively (if not unique, 
any solution is sufficient). These satisfy 

< I < (24) 


\ai\ = 

n 


( 21 ) 


we call it a uniform weight state and represent it by the 
bar graph in Fig. [3] In this case we would apply the 
diagonal Hamiltonian "H = Umax AT, where 


fix. 

2ir 


K = - 


0 

62 

2ir 


0 


Vo 0 


( 22 ) 


271- 


to the uniform state |unif) for a time tqc = 27r/i/5niax, 
which gives the desired target. 

Typically the target is not a uniform weight state, as 
represented in Fig. ID In this case we use the solution 


excluding the case where both < signs are identi¬ 
ties (which would violate the assumption that the 
target is nonuniform). Therefore 

2. Perform a phase shift C/diag = that brings 

the probability amplitudes and to the 
form = \ai^J and with 

|a*minl < |a*maxl- Apply SES Hamiltonian (P, 
where K is a diagonal matrix with = 

6'7mi„/37r and = (6»*„„/37r) - i, the other 

elements zero, and tqc = SttS/ g^nax- This phase shift 
is necessary to prepare the state for the next oper¬ 
ation. 

3. Implement a partial iSWAP t/gwap = from 

component Zmax to imin to bring the weight of Imin 
to the uniform value. 


|unif) — IFdiag (f^swapb(diag)M* * * (fAwapb^diag) 1 IV^) (^3) 

to the inverse problem of constructing the uniform state 
|unif) from the target [H, [T^. Each of the M steps 
in (na) consists of a pair of operations C/diag and C/gwap 
that move weight between a pair of components. After 


kiminP (25) 

n 

and leaving component I„iax with weight 

bimaxP |a*minP + I ^Imax ^ “ —■ ( 26 ) 
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Apply SES Hamiltonian ([IJ with = 

ATimax.imin = ^ all othor elements zero, and 
tqc = with V? given by 

I COS V? + |ai,„ax I sin = s/^. (27) 

There is always a solution with 0<(/3<7r/2. 

This completes the first step. 

If after the first step (17swapC4iiag)iIV') is a uniform 
weight state, it can be written in the form 

e“Ml) + e-"|2) + -.. + e“-|n) 

- 


is a uniform weight state, after which Wdiag is applied. 
The number of iterations required satisfies 

M<n-\. (31) 

This completes the solution to the inverse problem (1231) . 
We now use (l23l) to obtain 

IV') = (4agt^slap)l ■ • • (t^LgC^.lap)M I Uuif) , (32) 

which solves the general state-preparation problem in 
0(n) steps. Hermitian conjugations are implemented by 
changing the signs of the K matrices given above. The 
protocol given in (|32ll is, by itself, practical for small n. 

The complete state preparation operation can be sum¬ 
marized as 

|V') = C/|1), (33) 


and we apply the final operation Wdiag = e to 

produce (I^Ul) . Here we use SES Hamiltonian © with 

0 ... 0\ 

0 ... 0 

K= . , . (29) 

VO 0 ... 

and tqc = 27rli/gmax. If (I7swapt^diag)i IV') is not a uniform 
weight state, we again find the minimum and maximum 
weight components imin and imax, and follow the above 
protocol to generate (t7swapIVdiag)2(IV3wapIVdiag)i|V'). The 
procedure is repeated until 

(bVwapIViiag)M * * ‘ (bVwapIVdiag )2(C^swapt/diag)l|V') (30) 

I 


where 

U ^ {uL,Ul,,h ■ ■ ■ (C/dVagC^slap)M Wl^ 

(34) 

is the compiled unitary of the state-preparation algo¬ 
rithm. The three-step state preparation protocol uses 
the ABA decomposition to implement (IMl) . The total 
state preparation time, not including the |1) state ini¬ 
tialization time, is given in (HU). 

For example, suppose we wish to prepare the randomly 
chosen target 

IV') = 0.4829 11) -f (-0.5478-0.0575ii) |2)-h (0.1142-^0.23871) |3) 
-h (0.4095-f 0.24001) |4) -h (-0.3215-f 0.25451) |5), (35) 

in the n = 5 graph, where for convenience the first com¬ 
ponent has been chosen to be real. Following the state- 
preparation protocol leads to the compiled unitary 


U = 


/ 0.4829 

-0.5478- 0.05751 
0.1142-^ 0.23871 
0.4095 -f 0.24001 
V-0.3215-h 0.25451 


0.4499- 0.01581 
0.5855- 0.41531 
0.4664 -h 0.07001 
0.0841 - 0.12711 
0.1071-f 0.15771 


0.4499-0.01581 
0.1778-0.02491 
-0.7862- 0.25821 
0.1471 - 0.14921 
0.1071-^0.15771 


0.4478 - 0.01331 0.3984 -H 0.04501 \ 

-0.1305-^0.27031 -0.0855-^0.22731 
0.0910- 0.02841 0.1145 - 0.02221 

-0.7941-^0.18181 0.1471 -0.14921 
0.1071-^0.15801 0.1399- 0.83861 / 


(36) 


up to a phase factor. The first column of (1361) is the and 
target state. The ABA decomposition (Hip then leads to 




/-3.0826 

1.8972 

0.3983 

0.8753 

0.5934 \ 

/-1.1145 0.1981 0.3247 

-0.0776 -0.1888\ 


1.8972 

-3.7784 

0.5761 

0.3537 

0.5581 

0.1981 -2.6988 0.0219 

-0.2069 -0.0249 

B = 

0.3983 

0.5761 

-3.2370 

0.1664 

0.2327 

0.3247 0.0219 -1.9798 

-0.5623 0.1052 


0.8753 

0.3537 

0.1664 

-2.6191 

0.1488 

-0.0776 -0.2069 -0.5623 

-0.5291 -0.0747 

V 0.5934 

0.5581 

0.2327 

0.1488 

-4.6171/ 


V-0.1888 -0.0249 0.1052 -0.0747 -1.7104/ (38) 

(37) The associated K matrices and evolution times are de- 
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termined from the procedure given in Sec. IIII Al 


Ka = 


9a = 1.0848, 


/ 0.4604 

0.1826 

0.2993 

-0.0715 

0.1826 

-1 

0.0202 

-0.1907 

0.2993 

0.0202 -0.3373 

-0.5183 

-0.0715 

-0.1907 -0.5183 

1 

\-0.1741 

-0.0229 0.0970 

-0.0689 


.1741\ 


0.0970 

-0.0689 


(39) 


and 


Kb 


9b 


/0.2822 

1 

0.2100 

0.4614 

0.3128 \ 

1 

-0.0845 

0.3037 

0.1864 

0.2942 

0.2100 

0.3037 

0.2009 

0.0877 

0.1226 

0.4614 

0.1864 

0.0877 

0.5266 

0.0785 

Vo.3128 

0.2942 

0.1226 

0.0785 

-0.5266/ 


1.8972. (40) 


The total state preparation time, not counting the |1) 
state initialization, is given by m- This is about 13 ns 
for the target state (1551) in an SES chip with gmax/27r = 
50 MHz. 

Although state preparation is implemented in three 
steps for any n, the runtime does have a weak n- 
dependence, because 9a and 9b do. Averaged over ran¬ 
dom targets we find that 


26U -h fe ~ 4.0 X (41) 


Standard readout of an SES processor consists of the 
simultaneous measurement of each qubit in the diagonal 
basis. The SES condition means that a single qubit will 
be found in the state |1), with the remaining n — 1 qubits 
in |0). Let i be the qubit observed in it’s excited state. 
The probability of observing the excitation in qubit i is 
Pi = (i|/9|i). Therefore, if we have access to multiple 
copies of p we can repeat the readout N times to obtain 
estimates of the occupation probabilities pi with sam¬ 
pling errors no larger than {2y/N)~^. 

To compute (O), perform a (classical) spectral decom¬ 
position to a unitary V containing the eigenvectors of O 
as columns, and a real diagonal matrix D: O — VDVK 
Then we have 

(O) = TTipVDV'f) = Tt{pD), (43) 

where 

p' = V^'pV. (44) 

Therefore we can compute (O) by applying the unitary 
operator using the ABA decomposition, measuring 
the resulting occupation probabilities, which we denote 
by to indicate the application of V\ and then clas¬ 
sically evaluating the quantity 

n 

{0) = Y,DuPr''- (45) 

i=l 


For small n, either the linear-depth protocol (1551) or the 
three-step protocol based on (1551) can be used. However 
for large n, only the three-step protocol is practical. 


C. Computation of expectation values 

Finally, we show how to compute the expectation value 

(O) = Tr(pO) (42) 

of any nxn Hermitian observable O, by implementing 
the protocol of Reck et al. [3. Here p is any pure or 
mixed SES state provided as an input to the procedure. 


V. CONCLUSIONS 

In this work we have extended the SES method of 
Ref. Q to include a three-step implementation of arbi¬ 
trary nxn unitaries. The fast state preparation proto¬ 
col of Sec. IIVBI should be especially useful for practical 
quantum computing applications. 
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